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Dynamical Chern-Simons (DCS) modified gravity is an attractive, yet relatively unexplored, can- 
didate to an alternative theory of gravity. The DCS correction couples a dynamical scalar field 
to the gravitational field. In this framework, we analyze the perturbation formalism and stabil- 
' ity properties of spherically symmetric black holes. Assuming that no background scalar field is 

present, gravitational perturbations with polar and axial parities decouple. We find no effect of the 
Chern-Simons coupling on the polar sector, while axial perturbations couple to the Chern-Simons 
s ! ' scalar field. The axial sector can develop strong instabilities if the coupling parameter f3, associated 

to the dynamical coupling of the scalar field, is small enough; this yields a constraint on (3 which is 
. much stronger than the constraints previously known in the literature. 
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I. INTRODUCTION 



Interferometric and resonant-bar gravitational-wave detectors are now working at or near design sensitivity. It is 
expected that instruments such as LIGO or any of its advanced versions 1] will soon make the first direct detection 
of gravitational waves on Earth. The detection and measurement of gravitational waves from compact, massive 
astrophysical bodies opens a new window into the universe, and it also opens up the exciting possibility of testing 
Einstein's field equations in an unprecedented way. From measurements of the inspiralling phase of black hole or 
neutron star binaries.one can test General Relativity's prediction for the waveform and among others, bound the 
mass of the graviton [2|, y, |4[ . Gravitational- wave observations of the ringdown phase of the final black hole can lead 
(~-« ■ to tests of the no-hair theorem in General Relativity [H, HI]. Taken together, these measurements also allow for tests 
' of Hawking's area theorem for black holes [7J ■ 

What if there are corrections to the field equations, and do they leave a measurable imprint on gravitational waves? 
Unfortunately, this important question depends on several unknowns. It depends on the form of the corrections to 
the field equations, on the black hole solutions to these modified equations and on which specific process is generating 
the gravitational wave signal. One of such theories is known as Einstein-Dilatonic-Gauss-Bonnet, a simple example of 
one-loop corrected four-dimensional effective theory of the heterotic superstrings at low energies [||. The possibility 
of astrophysical tests of this theory were studied in Ref. Q . 

Another promising extension of general relativity is Chern-Simons (CS) gravity fiol. ITU IT3|. in which the Einstein- 
Hilbert action is modified by adding a parity-violating Chern-Simons term, which couples to gravity via a scalar 
field. This correction arises in many contexts. Such a term could help to explain several problems of cosmology, from 
inflation (as discussed by Weinberg [TH) to baryon asymmetry [3, EH, E6|. In most of the moduli space of string 
theory, a CS correction is required to preserve unitarity; furthermore, duality symmetries induce a CS term in all 
string theories with a Ramond-Ramond scalar [T3| • In loop quantum gravity, it is required to ensure gauge invariance 
of the Ashtekar variables [l8|] and it also arises naturally if the Barbero-Immirzi parameter is promoted to a field 
011. 

An interesting feature of CS modified gravity is that it has a characteristic observational signature, which could 
allow one to discriminate an effect of this theory from other phenomena. Indeed, the CS term violates parity, and thus 
it mainly affects the axial-parity component of the gravitational field. For instance, it yields amplitude birefringence 
of gravitational waves (the two polarizations travel with the same speed, but one is enhanced, while the other is 
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suppressed); on the other hand, the Schwarzschild solution is unaffected by CS modified gravity, and then the solar 
system tests of general relativity do not put strict bounds on the magnitude of this correction. Another possible 
signature of CS modified gravity may be found in the polarization of primordial gravitational waves [2ll |. 

Most of the literature on CS modified gravity refers to its non-dynamical formulation, in which the scalar field is a 
prescribed function. Usually, the so-called canonical prescription 9 oc t is chosen [22|, HHHij]- An explicit prescription 
for the scalar field is necessarily a priori, and furthermore it breaks gauge invariance (see the discussion in [25j]). For 
these reasons, we prefer to consider the dynamical Chern-Simons (DCS) modified gravity [26], where the scalar field 
is treated as a dynamical field. Note that CS gravity and DCS gravity are inequivalent and independent theories: 
although the CS gravity action can be obtained as a certain limit of the DCS gravity action, the solutions of CS 
gravity cannot be obtained from the solutions of DCS gravity [25]. We also remark that CS modified gravity, with 
the canonical prescription, is retrieved in the case of weak gravitational field, if one only keeps the leading order 
perturbations of Minkowski spacetime [27| ; therefore, CS gravity is appropriate to study, for instance, deviations from 
general relativity in solar system astrophysical processes, or in the motion of binary pulsars far from coalescence, but 
DCS gravity is required to describe deviation from general relativity in the strong field regime, like in black hole (BH) 
perturbation theory or in the coalescence of binary systems. 

In this paper we study, in the context of DCS gravity, gravitational perturbations of a spherically symmetric black 
hole. Perturbations of black holes are interesting for several reasons: since black holes populate the universe in large 
number, assessing their stability in a given theory is tantamount to testing the theory: if the black hole solution 
in that specific theory is unstable, they would not be seen. Second because a variety of processes taking place in 
the vicinities of black holes will be observed by gravitational-wave detectors. The two most important are extreme- 
mass-ratio-inspirals [28| . and quasi-normal ringing [f|. The first consists, for instance, on a small star orbiting around 
a supermassive black hole. Such a process can be modeled as a test particle inducing perturbations on a black 
hole background, and can be tackled with perturbation tools. The second is a universal signal: all or almost all 
events involving black holes produce a gravitational wave signal which at late times consists of a superposition of the 
characteristic modes of the black hole, the quasi- normal modes (QNMs). Such a signal will be seen by present or 
future gravitational wave detectors. 

BH perturbations have been studied, in the context of CS gravity, in Ref. [29], where it was found that polar- 
parity and axial-parity gravitational perturbations are coupled, and the equations do not allow for generic black 
hole oscillations. BH perturbations in DCS gravity were also briefly discussed in that work, in the framework of two- 
parameter perturbation theory: the two parameters were e, describing the magnitude of the gravitational perturbation, 
and t, describing the magnitude of the CS scalar field. The authors of [29| found that, if the scalar field has both 
a background (spherically symmetric) component of order 0(t) and an oscillating component of order 0(re), polar 
and axial perturbations are coupled, and the equations are extremely involved. Here, we consider the case in which 
the background scalar field vanishes, i.e. the scalar field is only generated by gravitational perturbations. Under 
this condition, we find that polar and axial gravitational perturbations decouple, and only axial parity gravitational 
perturbations are affected by the CS scalar field. Furthermore we find that, under this assumption, gravitational 
perturbations and the Chern-Simons scalar field are described by a simple set of equations. Numerical integration of 
these equations is not an easy task, due to the well-known asymptotic di verg ence which prevented for many years the 
numerical computation of quasi-normal modes of the Schwarzschild BH (30, |31|. 

We find some evidence that in the limit (3 — > oo, where (3 is the coefficient in front of the kinetic CS term of the action 
(see Eq. (jTJ) and discussion below), the BH does not admit QNMs. Furthermore we find that for < (3 M 4 < 2tt, 
there is at least one (strongly) unstable mode, and the Schwarzschild DCS solution is then unstable. This seems to 
impose a strong constraint on the theory, if it is to be compatible with the existence of astrophysical black holes, i.e., 
j3 > 1CU 2 km~ 4 . Previous bounds were much weaker: in [2 51 ] it was found that observational data from the extreme 
double pulsar system PSR J0737-3039 imply (3 > KT 15 km" 4 . 

Strictly speaking, our results apply only to non-rotating BHs, since rotating BH have a non-vanishing background 
scalar field |25l . [HJ and are more complex to handle. Astrophysical BHs are in general rotating (perhaps even rapidly 
rotating), but it is reasonable to expect that the range of [3 incompatible with the existence of non-rotating BH is 
also incompatible with (slowly) rotating BH. Therefore, we take the bound (3 > 10~ 2 km -4 as a (strong) indication 
of the values we may expect for the parameter (3. Finally, we also want to stress that the bound (3 > 1CU 2 km~ 4 does 
not at all rule out the possibility of an observational signature from DCS gravity. Indeed, we restrict the range of /3, 
but we do not restrict neither the range of a (i.e., with the a = 1 normalization, the amplitude of the scalar field), 
nor the time derivative of the scalar field amplitude. For instance, in (27| it was found that observational data from 
the double pulsar J0737-3039 imply a bound on i9, in the context of CS gravity. 

The paper is organized as follows. In Section[TT]we briefly review DCS modified gravity. In Section HTT1 we derive the 
equations for perturbations of a spherically symmetric black hole in DCS modified gravity. In Section [TV] we discuss 
the perturbation equations we have derived, finding some of their solutions. In Section [V] we draw our conclusions. 
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II. DCS GRAVITY 

Following the notation of [251 ] . the action of dynamical Chern-Simons modified gravity is 

S = k J d 4 x^R + | J d^xj—g-d *RR-^J d^x^g [g ab W a tfV 6 tf + V(0)] + S rnat . (1) 

Note that there are two parameters a, (3, but one of them can be eliminated by choosing the normalization of the 
scalar field. For instance, in [2t| the normalization of $ is chosen such that = 1, and a is the only parameter 
describing the coupling between the scalar field and the gravitational field. We normalize instead the scalar field by 
imposing a = 1, and use geometrical units c = G = 1 so that k — j^. Then, 

{S}=1 2 , [d] = l 2 , [/3]=r 4 . (2) 

Furthermore, we neglect V($), and consider the vacuum solutions (S ma t = 0); therefore, the equations of motion are 



R ab = -167TC.6 + 87r£0,„1?,6 (3) 

1 
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a-d = — — *RR (4) 



where 

C Qb = §, a e cde{a \I e R b) d + fi, dc *R d ( ab ) c (5) 
*RR = l -R abcd e ab ^ R cd ef . (6) 

Equations (J3J), ([3]) acquire a particularly simple form in the spherically symmetric case [HHHI, yielding (up to 0(d 2 )) 
the Schwarzschild solution; indeed, in the Schwarzschild spacetime *RR = and, assuming that the scalar field is also 
spherically symmetric § = i?(t, r), C ab — 0. Then, if we neglect the CS stress-energy tensor (quadratic in the scalar 
field), Eqns. ©, (01) are satisfied by the Schwarzschild solution. 

We mention that in [25|, [52] the solution for slowly rotating black holes has been found in DCS modified gravity; 
the corrections from the general relativistic solutions are of the order a//3, a 2 /(/3k), i.e., with our normalizations, of 
the order /3 _1 . Furthermore, an observational constraint for /3 _1 has been derived in [25| . from frame dragging effects 
in the extreme double pulsar system PSR J0737-3039 A/B: 

/T 1 < 10 15 km 4 . (7) 



III. PERTURBATIONS OF A SCHWARZSCHILD BACKGROUND 



We now consider perturbations of the spacetime geometry away from a Schwarzschild background. We set gffl to 
be the Schwarzschild metric and choose the Regge- Wheeler gauge for perturbations: 
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where Y are the scalar spherical harmonics and 
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(9) 



Here, (Ho, Hi, H-i, K) lm are (functions of r) describing the polar parity metric perturbations, (ho, h\) lm describe the 
axial parity metric perturbations. 

We now make one simplifying, self-consistent assumption, namely that the Chern-Simons scalar field is of the order 
of 0(h). In other words, the background scalar field the solution of the homogeneous equation n$(°> = 0, 
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is vanishing: = 0. The only scalar field present is induced by the perturbations. In this respect the case we 
are considering is different from that considered in Section VI of Ref. [29]. Under this assumption, the harmonic 
expansion of the scalar field i? is 



r 



(10) 



From here onwards, we will drop the lm superscripts. 



Replacing expansions ([8]), (fT0|) into Eqs. ([3]), ([4]) and neglecting terms quadratic in the metric perturbations or 
in the scalar field, we find that perturbations with different parities decouple: polar parity metric perturbations are 
unaffected by the scalar field. These will not be discussed any further here, since the stability properties and QNMs 
of these solutions are very well known (see [33j and references therein; for recent reviews see [g, H3|). On the other 
hand, the axial parity metric perturbations are coupled with the scalar field. Indeed, equations ((3]) yield 



E x = r 3 (-AM + l(l + l)r)h - rf (2ir 4 uj hi + 192tt MQ + ir 5 Luh\ - 96tt MrO' 
E 2 = -iur 3 (2h - iruj hi - rh' ) + r 2 f (I 2 + I - 2) hi - 96tt i MwQ = , 
E 3 ee ir 3 uj h + rf (+2Mhi + r 2 fh\) = . 

We note that these three equations are not all independent, it is easy to show that 

/r 4 



(E 2 /r 



E3+ d-m + 2)r Ei 
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Combining (JTTJ, $TZ§ and |T 
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we finally get 
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and defining the Regge- Wheeler master function Q(r) by 
hi = f- l rQ. 
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(11) 
(12) 

(13) 



(14) 

(15) 
(16) 

(17) 
(18) 



Equations (fT7|) and (| 18|) form a system of coupled second order differential equations for the perturbations Q , 
Q lm (with dimensions [Q lrn ] — 1°, [Q lm ] — I 3 ), from which one can completely characterize the axial parity metric 
perturbations and the scalar field. 



IV. INTEGRATION OF THE PERTURBATION EQUATIONS 

Despite their apparent simplicity, numerical integration of the perturbation equations (|17p , (|18| is not an easy task. 
We first note that, as r* — > ±oo (i.e. r — * r H = 2M, r — * +oo), they reduce to the simple wave equations 

,^ + - 2 ) e "" = ($ + ^)°"" = 0; (19) 

then, every solution (Q, <d) of the perturbation equations has the asymptotic form, as r* — > ±oo, 

V — A H,oo e ' A H,oo e > 

9 = Bf^e 1 ""-- + Sg i00 e- i "'-* , (20) 

where w, -A^J 1 ^ 11 , B ^^ are complex numbers. A QNM is a solution of the perturbation equation which satisfies the 
Sommerfeld boundary conditions 

A^ = B°^=X» =BZ = 0, (21) 
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i.e., no radiation outgoing from the horizon, no radiation ingoing from infinity. 

By fixing the normalization — 1, and defining ^4'^ — A € G , we have a unique solution of equations (| 1 7[) and 
(fTSf : thus we end with two complex conditions (A 1 ^ = B™ = 0) to satisfy, and two complex numbers to determine 
(Ao and lo). This naive counting of degrees of freedom tells us that, as in the case of the Schwarzschild BH, the 
dimension of the space of solutions is zero, allowing for a discrete set of complex values for uj: the quasinormal modes 
of the black hole. The difference from the Schwarzschild case is that here there are two coupled differential equations, 
instead of a single Schroedinger-like equation. 

We are not aware of any discussion in the literature regarding the QNMs of coupled systems of ODEs such as 
the one we have here. The main difficulty in integrating equations (fT"7f and (fT8|) is known: for stable spacetimes 
Im(u>) < and the correct behavior at (say) infinity is e 1 " 1 **, which is exponentially dominant over the unwanted 
e -\wr, _ j n case f ^ c Schwarzschild BH, many routes have been explored to circumvent the divergence problem 
discussed above; the main successful approaches are: a reformulation of the equation as a recurrence relation, which 
is then expressed as a continued fraction [35|: a WKB approximation of the potential 36] ; an analytic continuation to 
the complex plane of the radial coordinate [33] . It seems quite difficult to extend any of these approaches to a system 
of two coupled equations like (jTTJ) and (|18p . For instance, if one reformulates our perturbation equations in terms of 
recurrence series (see Appendix [XJ , one finds two coupled series, with ten terms at each order, which do not seem to 
be expressible in terms of continued fractions. 

Therefore, we were not able to perform a full, numerically accurate search for the QNMs of the coupled system 
(jTTJ) and (fT5|) . We studied only two limiting regimes: (3M 4 3> 1 and /3M 4 < 1. Before proceeding to the analysis of 
these two cases, we mention that, since Q, © exponentially diverge in the r* — > ±oo limits, it is not at all obvious 
that equations (jTTJ) , (fT5|) reduce to the form (|19p and then that their solutions, in these limits, have the form (|20p . 
as we have assumed at the beginning of this section. Actually, Eq. (f2"0"l) can be proved without assuming (fT5|) . as we 
show in Appendix [Bl 



A. Very large ft 



In the case /3Af 4 > 1 Eqns. JI7J, (jTgJ) reduce to 



——D ln 



1(1 + 1) 6M 



Qlm = _ f 9_^M Qlr , 



e !m = o. 



(22) 
(23) 



Thus, in this limit, the eigenvalue problem for © is homogeneous, and exactly equivalent to the Schwarzschild case. 
The solutions of (|23[) . u>Q NM , 8^™, are well-known and can be computed accurately using the continued fraction 
approach 0, [35[ . Then, we are left with equation (j22|) , which can be considered as a single differential equation with 
source. One can build the solution to the inhomogeneous problem (1221) by using the Green function method: 



f 

W 



QnSdr + Q H / Q^Sdr 



2M 



(24) 



where Qh, Qoo are two independent solutions of the homogeneous version of (|22| (evaluated with the w 



eigenvalue of the scalar equation), such that Q c 
W is the wronskian between Qh, Qoo, and S : 



at infinity and Qh ~ e 



-e QJVM 967riwM/r 



close to the horizon; the quantity 



Solution (|24|) can be seen to satisfy the required boundary conditions and indeed is the only solution satisfying such 
boundary conditions. Unfortunately, the integrals in (|24D are not well defined: the wavefunctions Qh, Qoo and S 
diverge exponentially at both the horizon and spatial infinity. This seems to indicate that this solution does not exist, 
i.e. that DCS Schwarzschild BH do not admit QNM in which both the gravitational perturbations and the scalar 
field are excited. Only trivial QNM solutions, in which O = and the gravitational field oscillates at the frequencies 
of ordinary Schwarzschild QNM, seem to be allowed in this limit. 



B. Instability regime 



The opposite regime concerns very small [3. As might be anticipated from the previous discussion, the behavior of 
the solutions of the coupled system are intimately connected with the behavior of the homogeneous problem. Now, 
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in the case of small /?, the equation for the scalar field is 

with V = f (1(1 + l)/r 2 (l - 5767rM 2 /?~ 1 /r 6 ) + 2M/r 3 ). A study of the homogeneous equation (i.e., we set Q = 
in the above), shows that unstable modes are possible since the effective potential is negative in some regions. Since 
the potential is bounded in — oo < r* < +oo and it vanishes at both extrema, a sufficient condition for existence of a 
bound state of negative energy is [H, Hi| 



7^.)*. = 7 ^' 2 + 2 ' + 2 ^ t - 18 "' + 1K <°. M 

This condition is satisfied whenever 

18Z(Z+ 1)7T 

P 7(2/ 2 + 21 + 1)M 4 ' 1 J 

which yields (3 < 108tt/(91M 4 ), 216tt/(175M 4 ) for I = 2, 3 respectively. An example of an unstable I = 2 mode of the 
homogeneous equation is shown in Figure [TJ for (3=1. This does not prove the existence of unstable modes for the 
coupled system, but it is a strong argument in its favor. 

Numerical integration of the coupled system (TIT)) . (fl~8"|) shows that unstable modes, indeed, do exist. In our search, 
we set Btf = 1, and look for solutions with lo purely imaginary; in this case, the boundary conditions (|2ip imply that 
the functions Q, vanish at r — > oo. 



0,04 




FIG. 1: Left: Scalar wavefunction O for the unstable mode corresponding to M 4 f3 = 1 and Mui — 0.51197595i, assuming no 
gravitational coupling, i.e., this is an unstable mode of the homogeneous version of eq. (|25[) . Right: Scalar wavefunction Q and 
Regge- Wheeler function Q for the unstable mode of the full coupled system, corresponding to M 4 (3 = 1 and Muj = 0.4280636i, 
AS = -0.1069141. 



A typical unstable mode of the coupled system is shown in FigureJTJ for (3 = 1. The overall qualitative and quantitative 
behavior of the scalar field is in agreement with the analysis of the homogeneous scalar wave equation. The 
instability timescale increases (i.e. Imw decreases) when f3 increases. For instance, Muj ~ 0.428i, 0.263i for M i /3 = 1,2 
respectively. This was to be expected from the analysis of the homogeneous scalar equation. We were not able to find 
unstable modes for f3M A larger than ~ 2n. This is again in agreement with the heuristic arguments above, though a 
more systematic approach to investigate this issue is necessary. 

We have then shown that spherically symmetric BH in DCS modified gravity admit unstable modes for (3M A < 27r; 
therefore, since astrophysical black holes exist with masses M astro ^ 3Af Q ~ 5 km, we have the constraint 



(28) 
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much stronger than (7}. The present analysis does not directly apply to astrophysical BHs, since these are rotating 
in general. Rotating black holes in this theory require a non-vanishing background scalar field pBl |32| and therefore 
our entire analysis needs to be modified. On the other hand, it is hard to imagine that rotation can drastically alter 
the instability regime. We thus believe that the bound l|28p should be considered as a strong indication on the values 
that we may expect for (3 in the general case. 

We remark that the normalization choice a = 1 does not really affect our results. If we do not fix the parameter a, 
some terms in equations (IT71) . (TT5)) get multiplied by a or a 2 , and the bound (|28p becomes 

£>^>10-W. (29) 
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Note that the quantity a 2 /{(3k) is what is called £ in [25|, where the bound Q is derived. 

V. CONCLUSIONS 

We have studied perturbations of black holes in the context of dynamical Chern-Simons theory. Under the as- 
sumption that the scalar field is of the order of the gravitational perturbation, we found that axial and polar parity 
gravitational perturbations are decoupled, and only axial parity perturbations are coupled with the scalar field. The 
equations describing these perturbations are fairly simple, but their numerical integration is tricky: our attempts to 
generalize to DCS gravity either the recurrence relation approach of [HI, |4(j, or the Riccati equation approach of 
[30(, have been unsuccessful. Still, we have been able to extract some interesting information from the perturbation 
equations. 

We found that black holes in this theory can develop strong instabilities. The constraint on the coupling parameters 
necessary to avoid the instability is stringent. We obtain (3M A > 2ir for the spacetime to be stable, which translates 
to j3 > I CP 2 km~ 4 . Thus, the observation of stellar-mass black holes imposes a constraint on the coupling parameter 



which is about 10 13 times more stringent than bounds from binary pulsar dragging effects (Yunes and Pretorius 25 1 
obtain (3 > f CP 15 km~ 4 from the pulsar PSR J0737-3039 A/B). This bound has been derived under the assumption 
that the BH is non-rotating, but it should provide a reliable estimate on the range of allowed in the general case. 

Furthermore, we found some evidence that in the limit (3 — > oo, there are no non-trivial QNM solutions: the only 
allowed solutions in this limit seem to be the ordinary Schwarzschild gravitational QNM, with a vanishing scalar field. 

Much more remains to be done. The coupled system of equations (fT7|) and (TT5]) was not solved in general in the 
present paper, nor were the exact limits of instability investigated. It would be extremely interesting to study the 
spectra of this system, either by analytical approximations or by numerical means. A promising approach would 
consist in doing wave scattering using time-evolution methods [4l|, and reading off the quasinormal modes directly 
from the frequency and decay time of these perturbations. 

The formalism developed here can be used to study point particles in the background geometry and to compute 
the gravitational wave signal generated by extreme-mass-ratio inspirals. Generalization to rotating black holes is a 
major (rotating black hole solutions are only partially understood [12j|). but fundamental task to accomplish. 
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APPENDIX A: RECURRENCE RELATIONS 

Here we try to generalize the continued fraction method [H, to solve Eqns (JTTJ) , (|I8p . 
The functions Q, O with the correct asymptotic behavior can be written as 

Q = (— l) —) e ^-2M )q> 



\2M J \2MJ 

J- -A iui(r — 2M) \ 



6 = — -1 — e'^-^E, (Al) 

\2M J V 2M J ' v ; 
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where E can be expressed as 

oo oc 

$ = any 71 , E = Ky n (A2) 

n=0 n=0 

in terms of the dimensionless variable 

y = 1 . A3 

r 

Note that a n< o = 6 n <o — 0; we set the overall normalization by imposing ao = 1, and the (complex) constant &o is 
an unknown of the problem. 

Replacing these expressions in the perturbation equations (|17|) . {T5J we find the recurrence relations 

+ A (& n _ 2 CT* + . . . + fr n _7(j£) 

i/3„ + + p + . . . + a„_ 7 cr^) 



(A4) 



where 



c*n — n ( n — 4iMw) 

/3„ = -2n 2 + (16iA/w + 2)n - 8iAfw + 32A/V - Z(Z + 1) - 1 
j n = n 2 - (8iMw + 2)n + 8iMw - 16AZ V + 1 
9ttZ(Z + 1) 



7n = 7n - 4 - 
5„ 



A 



/3Af 4 
54tt/(7 + 1) 

/3M 4 

135ttZ(Z + 1) 

/3M 4 
37riw 



M 4 

(Z + 2)! 3iw 
(Z - 2)! 16/3 A/ 4 

(1,-5,10,-10,5,-1). (A5) 



The two complex recurrence relations (|A4[) . depending on the two unknown complex quantities (to, bo), in principle 
converge only for a discrete set (w, &o)j, corresponding to the QNM; it is not clear how to numerically implement this 
condition, since it seems not possible to express the recurrence relations (|A4I) in terms of continued fractions. 



APPENDIX B: ASYMPTOTIC BEHAVIOR OF THE PERTURBATION EQUATIONS 



We first note that, as r* — > ±oo (i.e. r — > rjj = 2M, r — > +oo), the perturbation equations (jTTJ) and (fT5| reduce to 
the simple coupled wave equations 



Changing scalar wavefunction 



+ w 2 J e 

+ lu 2 ) Q 



<ir 2 
dr 2 



f (l + 2)\6iM 
'r*(l - 2)! ujf3 ' 
f 



(* + 2)! 
4u;V (Z - 2)!tt/? 



0" 



(Bl) 
(B2) 

(B3) 
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we get 



^ + - 2 )e- = -^Q, (B4) 



with T = 2UMsJiTfii. Finally, denning 4<+ = Q + 9 ncw and = Q - 6 ncw , one gets the equivalent system 



fl=5JT 

+ co 2 *+ = -^*+, (B6) 



dr 2 

^ + J) 9- = tl*-. (B7) 



Each of the wavefunctions 4 ,± has a simple harmonic behavior at the boundaries and therefore ([^0]) follows 
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